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Abstract 

The road coloring theorem states that every aperiodic directed 
graph with constant out-degree has a synchronized coloring. This the- 
orem had been conjectured during many years as the road coloring 
problem before being settled by A. Trahtman. 

Trahtman's proof leads to an algorithm that finds a synchronized 
labeling with a cubic worst-case time complexity. We show a variant 
of his construction with a worst-case complexity which is quadratic in 
time and linear in space. We also extend the road coloring theorem to 
the periodic case. 



1 Introduction 

Imagine a map with roads which are colored in such a way that a fixed 
sequence of colors, called a homing sequence, leads the traveler to a fixed 
place whatever be the starting point. Such a coloring of the roads is called 
synchronized and finding a synchronized coloring is called the road coloring 
problem. In terms of graphs, it consists in finding a synchronized labeling 
in a directed graph. 

The road coloring theorem states that every aperiodic directed graph 
with constant out-degree has a synchronized coloring (an aperiodic graph 
is strongly connected and the gcd of its cycles is equal to 1). It has been 
conjectured under the name of the road coloring problem by Adler, Goodwin, 
and Weiss [2], and solved for many particular types of automata (see for 
instance [2], [8], [7], [lO]). Trahtman settled the conjecture in [TT]. 

The problem is not only a puzzle but has many applications in various 
areas like coding or design of computational systems. That systems are of- 
ten modeled by finite-state automata (i.e. graphs with labels). Due to some 
noise, the system may take a wrong transition. This noise may for instance 
result from the physical properties of sensors, from unreliability of compu- 
tational hardware, or from insufficient speed of the computer with respect 
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to the arrival rate of input symbols. It turns out that the asymptotic behav- 
ior of synchronized automata is better than the behavior of unsynchronized 
ones (see f^). Synchronized automata are thus less sensitive to the effect of 
noise. 

In the domain of coding, automata with outputs {i.e. transducers) can 
be used as encoders or decoders. When they are synchronized, the behavior 
of the coder or of the decoder is improved in the presence of noise or er- 
rors. For instance, the well known Huffman compression scheme leads to a 
synchronized decoder provided the lengths of the codewords of the Huffman 
code are mutually prime. It is also a consequence of the road coloring the- 
orem that coding schemes for constrained channels can have sliding block 
decoders and synchronized encoders (see [l] and [9]). 

Trahtman's proof is constructive and leads to an algorithm that finds 
a synchronized labeling with a cubic worst-case time complexity [12]. The 
algorithm consists in a sequence of flips of edges with the same origin so 
that the resulting automaton is synchronized. One first searches a sequence 
of flips leading to an automaton with a stable pair of states {s,t) (i.e. with 
good synchronizing properties). One then quotients the automaton by the 
congruence generated by (s, t) and iterates the process on this smaller au- 
tomaton. Trahtman's method for finding the flips leading to a stable pair has 
a worst-case quadratic time complexity, which makes his algorithm cubic. 

In this paper, we design a worst-case linear time algorithm for finding 
a sequence of flips until the automaton has a stable pair. This makes the 
algorithm quadratic in time and linear in space. The sequence of flips is 
obtained by fixing a color, say red, and by considering the red cycles formed 
with red edges, taking into account the positions of the root of red trees 
attached to each cycle. The prize to pay for decreasing the time complexity 
is some more complication in the choice of the flips. 

The article is organized as follows. In Section [21 we give some definitions 
to formulate the problem in terms of finite automata instead of graphs. 
In Section [3] we describe Trahtman's algorithm and our variant. We give 
both an informal description of the algorithm with pictures illustrating the 
constructions, and a pseudocode. The complexity is analyzed in Section [5l 

2 The road coloring theorem 

In order to formulate the road coloring problem we introduce the notion 
of automaton. A (finite) automaton A over some (finite) alphabet A is 
composed of a finite set Q of states and a finite set E of edges which are 
triples (p, a, q) where p, q are states and a is a symbol from A called the label 
of the edge. 

An automaton is deterministic if, for each state p and each letter a, 
there is at most one edge starting in p and labeled with a. It is complete 
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deterministic if, for each state p and each letter a, there is exactly one edge 
starting in p and labeled with a. This implies that, for each state p and each 
word w, there is exactly one path starting in p and labeled with w. If this 
unique path ends in a state q, we denote hy p ■ w the state q. 

An automaton is irreducible if its underlying graph is strongly connected. 
It is aperiodic if it is irreducible and if the gcd of the cycles of its graph is 
equal to 10. 

A synchronizing word of a deterministic complete automaton is a word 
w such that for any states p, q, one has p ■ w = q ■ w. A synchronizing word 
is also called a reset sequence or a magic sequence, or also a homing word. 
An automaton which has a synchronizing word is called synchronized (see 
an example on the right part of Figure [T]) . 

Two automata which have isomorphic underlying graphs are called equiv- 
alent. Hence two equivalent automata differ only by the labeling of their 
edges. In the sequel, we shall consider only complete deterministic automata. 

The road coloring theorem can be stated as follows. 

Theorem 1 (A. Trahtman Any aperiodic complete deterministic au- 

tomaton is equivalent to a synchronized one. 




Figure 1: Two complete aperiodic deterministic automata labeled on ^ = 
{a,b}. A thick red plain edge is an edge labeled by a while a thin blue 
dashed edge is an edge labeled by b. The automaton on the left is not 
synchronized. The one on the right is synchronized; for instance, the word 
aaa is a synchronizing word. The two automata are equivalent. 



3 An algorithm to find a synchronized coloring 

Trahtman's proof of Theorem[T]is constructive and gives an algorithm to find 
a labeling that makes the automaton synchronized provided it is aperiodic. 

In the sequel A denotes an n-state deterministic and complete automaton 
over an alphabet A. We fix a particular letter a (z A. Edges labeled by a 
are also called red edges. The other ones are called b- edges. 

A synchronizable pair of states in an automaton is a pair of states {p, q) 
such that there is a word w with p ■ w = q ■ w. It is well known that an 
automaton is synchronized if all its pairs of states are synchronizable. 

^Note that this notion, which is usual for graphs, is not the notion of aperiodic automata 
used elsewhere and which refers to the period of words labeling the cycles (see e.g. [6]). 
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A stable pair of states in an automaton is a pair of states (p, q) such 
that, for any word u, (p- u,q ■ u) is a synchronizable pair. In a synchronized 
automaton, any pair of states is stable. Note that if (p, q) is a stable pair, 
then for any word n, {p-u,q ■ u) also is a stable pair, hence the terminology. 
Note also that, if (p, q) and (g, r) are stable pairs, then (p, r) also is a stable 
pair. It follows that the relation defined on the set of states by p = g if 
(p, q) is a stable pair, is an equivalence relation. It is actually a congruence 
{i.e. p ■ u = q ■ u whenever p = q) called the stable pair congruence. More 
generally, a congruence such that any pair of states in a same class is stable, 
is called a stable pair congruence. The congruence generated by a stable pair 
(s, t) is the least congruence such that s and t belong to a same class. It is 
a stable pair congruence. 

The computation of the congruence generated by (s, t) can be performed 
by using Union/Find functions computing respectively the union of two 
classes and the leader of the class of a state. After merging two classes 
whose leaders are p and q respectively, we need to merge the classes oi p ■ £ 
and q ■ i for any i £ A. A pseudocode for merging classed in given in 
Procedure Merge below. 

Merge( stable pair (s,i)) 

1 x <— Find(s) 

2 2/ <- FiND(t) 

3 ii X ^ y 

4 then Union(x, y) 

5 tor £eA 

6 doMERGE{x ■ i,y ■ i) 

If ^ = {Q,E) is an automaton, the quotient of ^ by a stable pair 
congruence is the B whose states are the classes of Q under the congruence. 
The edges of B are the triples (p, £, q) where {p, q) is an edge of A. The 
automaton B is complete deterministic when A is complete deterministic. 
The automaton B is irreducible (resp. aperiodic) when A is irreducible (resp. 
aperiodic). 

Lemma 2 (Culik et al. [1]). // the quotient B of an automaton A by a 
stable pair congruence is equivalent to a synchronized automaton B' , then A 
is equivalent to a synchronized automaton A' . 

Proof. Let B' be a synchronized labeling of B. The equivalent automaton 
A' is defined as follows. If there is an edge {p, £, q) in A but no edge {p, £, q) 
in B' , then there is an edge {p, c, q) in B' with I ^ c. We flip the labels of the 
two edges labeled i and c going out of p in A. This definition is consistent 
with the definition of the edges of B. After all flips, B' is still a quotient 
oiA!. 

Let us show that A' is synchronized. Let tt; be a synchronizing word 
of B' and r the state ending any path labeled by w in B' . Let p, q be two 
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states of A! . Then p ■ w and q ■ w belong to the congruence class of r. Hence 
(p ■ w, q ■ w) is a stable pair of A' . Therefore (p, q) is a synchronizable pair 
of A! . All pairs of A! being synchronizable, A! is synchronized. □ 

Trahtman's algorithm for finding a synchronized coloring of an aperiodic 
automaton A consists in finding an equivalent automaton A! of A that has 
a stable pair (s, i), then in recursively finding a synchronized coloring B' for 
the quotient automaton B by the congruence generated by (s, t), and finally 
in lifting up this coloring into the initial automaton as follows. If there is an 
edge ip,l,q) in A but no edge {p,i,q) in B' , then there is an edge (p,c,q) 
in B' with t ^ c. Then we fiip the labels of the two edges labeled i and c 
going out of p in A. The algorithm is described in the following pseudocode, 
where the procedure FindStablePair, which finds a stable pair of states in 
any aperiodic automaton of size at least 2, is described in the next section. 

FindColoring( aperiodic automaton A) 

1 A 

2 while A has at least two states 

3 do A, (s, t) <— FindStablePair(^) 

4 lift the coloring of A up to .4,0 

5 A'^ the quotient of A by the congruence generated par (s, i) 

6 return Aq 

The termination of the algorithm is guaranteed by the fact that the 
number of states of the quotient automaton of A is strictly less than the 
number of states of A. 

4 Finding a stable pair 

In this section, we consider an aperiodic complete deterministic automaton 
A over the alphabet A. We design an algorithm for finding an equivalent 
automaton with a stable pair, which has a linear-time complexity. 

In order to describe the algorithm, we give some definitions and nota- 
tions. 

Let H be the subgraph of the graph of A made of the red edges. The 
graph 7?. is a disjoint union of connected component called clusters. Since 

each state has exactly one outgoing edge in TZ, each cluster contains a unique 
(red) cycle with trees attached to the cycle at their root. If r is the root of 
such a tree, its children are the states p such that p is not on the a red cycle 
and {p, a, r) is an edge. If p, q belong to a same tree, p is an ancestor of q 
(or g is a descendant of p) in the tree if there is a red path from q to p. 

For each state p of the cluster, we define the level of p as the distance 
between p and the root of the tree containing p. If p belongs to the cycle, 
its level is null. For each state we denote by rooi[p] the root of the tree 
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containing p. A maximal state is a state of maximal level. A maximal 
cluster is a cluster containing a maximal state. A maximal tree is a tree 
containing at least one maximal state and rooted at a state of level 0. The 
level of an automaton is the level of a maximal state. A maximal root is the 
root of a maximal tree. 

Let C be the red cycle of a cluster and k its length. We number its states 
from to A; — 1 along the cycle. If p belongs to the cycle, we denote this 
number by num\p\. 

The algorithm relies on the following key lemma due to Trahtman 
It uses the notion of minimal images in an automaton. An image in an 
irreducible automaton A = {Q, E) is a set of states I = Q • w, where w is 
a word and Q-w = {q-w\q(z Q}. A minimal image in an automaton is 
an image which is a minimal element of the set of images for set inclusion. 
In an irreducible automaton two minimal images have the same cardinality 
which is called the minimal rank of A. 

Lemma 3 (Trahtman [H]). Let A be an irreducible complete deterministic 
automaton with a positive level. If all maximal states in A belong to the 
same tree, then A has a stable pair. 

Proof. Since A is irreducible, there is a minimal image / containing a max- 
imal state p. Let £ > the level of p {i.e. the distance between p and the 
root r of its tree). Let us assume that there is a state q ^ p G I of level £. 
Then the cardinal of / • a' is strictly less than the cardinal of I, which con- 
tradicts the minimality of /. Thus all states but p m I have a level strictly 
less than i. 

Let m be a common multiple of the lengths of all red cycles. Let C be 
the red cycle containing r. Let sq be the predecessor of r in C and si the 
child of r containing p in its subtree. Let J = I ■ a''^^ and K = J ■ a"^. 
Since the level of all states of / but p is less than I — 1, the set J is equal to 
{si} U R, where R is a set of states belonging to some red cycle. Since for 
any state q in a red cycle, q ■ a"^ = q, we get K = {sq} U R. 

Let whe a word of minimal rank (i.e. Q-w is a minimal image). For any 
word V, the minimal images J ■ vw and K ■ vw have the same cardinal equal 
to the cardinal of /. This implies that sq ■ vw = si ■ vw. As a consequence 
(so ■ v,si ■ v) is synchronizable and thus (sq, si) is a stable pair. □ 

We say that an automaton satisfies Condition (*) if all maximal states 
belong to the same tree. Note that Condition (*) implies that all maximal 
states belong to a same cluster. 

A set of edges going out of a state p is called a bunch if these edges all 
end in a same state q. Note that if a state q has two incoming bunches from 
p,p', then {p,p') is a stable pair. 

We describe below sequences of flips of edges that make the resulting 
equivalent automaton have a stable pair. 
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4.1 The case of null maximal level 

We first assume that the level of the automaton is i = 0. The subgraph 
TZ of red edges is a disjoint union of cycles. If the set of outgoing edges of 
each state is a bunch, then there is only one red cycle and the automaton is 
not aperiodic. Let p he a state whose set of outgoing edges is not a bunch. 
There exists h ^ a and q ^ r such that (p, a, q) and (p, 6, r) are edges. We 
flip these edges. This gives an automaton A which satisfies Condition (*). 
Let s be the state which is the predecessor of r on its red cycle. It follows 
from the proof of Lemma [3] that the pair (p, s) is stable. 

This case is described in the pseudocode FlipEdgesLevelZero where 
previous[q\ is the predecessor of q on its red cycle. The function FlipEdges- 
LevelZero(^) returns an automaton equivalent to A together with a stable 
pair. 

FLipEDGESLEVELZERO(automaton A) 

We assume that the maximal level of the states is 

1 for each state p on a red cycle C 

2 do 

3 if the set of outgoing edges of p is not a bunch 

4 then let {p, a, q) and (p, b, r) be edges with b a and q r 

5 flip these edges 

6 s ^ previous [q] 

7 return A, {p, s) 

8 return error(^ is not aperiodic) 

4.2 The case of non-null maximal level 
4.2.1 Preliminary flips 

Let us now assume that the level of the automaton is ^ > 0. We assume that 
the level of the states have been computed with a depth-first exploration of 
the subgraph TZ of red edges in a linear time. 

In a preliminary step, we perform a sequence of fiips of A which leads 
to an equivalent automaton satisfying the following Condition (**): for any 
state of maximal level p, there is no b-edge {t,b,p) where t is an ancestor of 
p distinct from the root in the tree containing p. 

In order to achieve this goal, we compute for each maximal tree T a set 
S of states t such that there is a 5-edge {t, b,p) with t ancestor of p inT, p 
of maximal level, t distinct from the root of T, and there is no other state 
t' ancestor of t satisfying the same condition. 

The sets 5 are computed in linear time by performing a depth-first ex- 
ploration of each tree. For each tree, we start at its root following the red 
edges in the backward direction and examining the 6-edges {t,b,p) in the 
forward direction. We assume that a function Exploration returns all 
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states t in 5 together with a 6-edge (t, b,p), p maximal, satisfying the above 
conditions. 

Wc then flip (t,b,p) and the red edge going out of t for each t € 5" as 
long as there still remains at least one maximal tree in the automaton. 

This sequence of flips is described in the function ScanTrees below. In 
the pseudocode, the value of previous[r,p] is the stable pair (so,si), where 
So is the predecessor of r on its red cycle, and si is the child of r which is 
an ancestor of p. 

Note that each such flip creates a new red cycle and hence a new cluster. 
The level of each state belonging to each new cluster decreases strictly. In 
particular, the level of p and t becomes 0. One can know whether there 
remains maximal states after each flip without recomputing all levels: with 
a depth-first search exploration of the trees, one computes for each state 
t S the number of maximal states descendant of t, and we update the 
whole number of maximal states after each flip. This can be done in a linear 
time complexity for the whole sequence of flips. 

The procedure either stops when we have flipped two outgoing edges of 
all states in S (and there still remains states of maximal level £) , or when we 
have flipped two outgoing edges of all states in S but one. In the first case, for 
any state of maximal level p, there is no 6-edge (t, b,p) where t is an ancestor 
of p distinct from the root of the tree of p {i.e. Condition(**) is fulfilled). 
In the second case, the resulting automaton A has a single maximal tree. 
Hence it satisfies Condition (*). Thus the function ScanTrees(^) returns 
either an equivalent automaton together with a stable pair of states, or just 
an equivalent automaton satisfying Condition (**). 

ScanTrees ( automaton A) 

We assume that the maximal level of the states is positive 

1 M •*— the number of states of maximal level 

2 5 <— Exploration (^) 

(to each state t G S is associated a 6-edge {t, b,p)) 

3 for each state t € S 

4 do n{t) ^ the number of maximal states descendant of t 

5 for each state t € S 

6 do 

7 if n{t) < M 

8 then flip {t, b,p) and the red edge going out of t 

9 M ^M- n{t) 

10 S^S-{t} 

11 if 5 = 

12 then return A 

13 else r ^ the root of the unique maximal tree of A 

14 return A, previous[r,p] 

The preliminary treatment is done only one time at the beginning and 
it has a linear time complexity. 
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4.2.2 Main treatment 

We now assume that the automaton satisfies Condition (**). We compute 
the level of the states with a depth-first exploration of the subgraph TZ of 
red edges (following the edges backward) in a linear time. For each node, 
one can access the root of his tree in constant time. 

Let C be a red cycle containing a maximal tree T rooted at r. Let sq 
be the predecessor of r in the red cycle. We denote by si,S2-, ■ ■ iSp the 
children of r which are ancestors of a maximal state, and by pi a maximal 
state descendant of Sj for 1 < i < p. We denote by tj a predecessor of 
Pi by a 6-edge {ti,bi,pi). We describe a procedure FlipEdges(^, r) which 
is a sequence of flips of edges such that the automaton obtained after the 
sequence of flips has a unique maximal tree. Note that the levels will not be 
recomputed after each flip (which would increase the time complexity too 
much). Each update will be explicitly mentioned. 

We consider several cases depending on the position of the states ti in 
A satisfying Condition (**). If there is a state ti which is not in the same 
cluster as r or if level[ti] > 0, then an automaton equivalent to A satisfying 
Condition (*) is obtained by flipping the edge {ti,bi,pi) and the red edge 
going out of ti. Indeed, one may easily check that, after the flip, all states 
of maximal level belong to the same tree as pi. 

If there is a state ti that belongs to the cycle C but which is not contained 
in the interval [r^ • . r[, where r, r2, . . is the ordered list of maximal roots, 
the same conclusion holds by flipping the edge {ti,bi,pi) and the red edge 
going out of ti , as is shown in Figure [2j 

Let us now assume that all states ti belong to the cycle C and ti G [r^ . . r[ 
for 1 < i < p. 

We first consider the case where there are at least two distinct states ti 
(see Figure E]) . Without loss of generality, we can for instance assume that 
ti < t2 in [rfc . . r[. We flip the edge (ii, bi,pi) and the red edge going out of 
ti. We denote by T' the new tree rooted at r. If the height of T' is greater 
than i, the automaton satisfles Condition (*) (see the right part of Figure[3|). 
Otherwise the height of T' is at most i (see the left part of Figure H]) . In 
that case, we also flip the edge {t2, ^2)P2) and the red edge going out of t2. 
The new equivalent automaton satisfles Condition (*) (see the right part 
of Figure H]) . The computation of the size of T' is detailed in the complexity 
issue. 

We now consider the case where there all states ti are equal to a same 
state t. This includes the case where the number p of children of r having a 
maximal descendant, is equal to 1. 

We flrst treat the case of p > 1 and all states ti equal to t. Let Tq be the 
tree rooted at r obtained if we flip {t,bi,pi) and the red edge going out of t, 
and if we keep only r and the subtree rooted at the child sq. The states of 
the tree Tq rooted at r are represented in salmon in the left part of Figure [5l 
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This step again needs a computation of the height of Tq explained in the 
complexity issue. If the height of Tq is greater than the height of T, we do 
the flip and the equivalent automaton satisfies Condition (*). If the height 
of To is less than the height of T, we also do the flip of {t,bi,pi) and the 
red edge going out of t. We then call again the procedure FlipEdges(^, r) 
with this new red cycle. This time, the height of the new tree Tq, denoted 
Tq is equal to the height of T. Hence this call is done at most one time for 
a given maximal root r. Finally, we consider the case where the height of 
T and Tq are equal (see the left part of Figure E]) • If the set of outgoing 
edges of sq is a bunch and there is a state Sj, for 1 < i < p, whose set of 
outgoing edges also is a bunch, we get a trivial stable pair (sq, s,). If the set 
of outgoing edges of sq is a bunch but none of the sets of outgoing edges of 
Sj for 1 < i < p is a bunch (as in the left part of Figure [5|), we flip (t, 
and the red edge going out of t. The height of the new tree Tq (obtained if 
we flip (i, 62,^2) and the red edge going out of t and keep only r and the 
subtree rooted at the child si) has the same height as T. We then call again 
the procedure FlipEdges(^, r) with this new red cycle. This time, the 
height of the new tree Tq is equal to the height of T and the set of outgoing 
edges of the predecessor of r on the cycle is not a bunch. This call is thus 
performed at most one time. 

If the set of outgoing edges of sq is not a bunch, let (sq, Qo) be a 6-edge 
going out of So with qq ^ r. If go does not belong to T, we get an equivalent 
automaton satisfying Condition (*) by flipping (sqj^OjQ'o) and the red edge 
going out of So. If qQ belongs to T, we flip (sq, ^O; 9o) and the red edge going 
out of So. We also flip (t, 61, pi) and the red edge going out of t if qQ is not a 
descendant of si, or (t, b2,P2) and the red edge going out of t in the opposite 
case. Note that sq ^ t since the height of Tq is equal to the non-null height 
of T. We get an equivalent automaton satisfying Condition (*) (see the 
right part of Figure E]) . 

We now treat the case of p = 1. As before we denote by Tq the tree 
rooted at r obtained if we flip (t, 61,^*1) and the red edge going out of t and 
keep only r and the subtree rooted at the child so. If the height of Tq is 
greater than the height of T, we do the flip and the equivalent automaton 
satisfies Condition (*). If the height of Tq is less than the height of T, we 
do not flip the edge {t,bi,pi) and the red edge going out of t, and return 
the automaton together with the edge {t,bi,pi). We now come to the case 
where the height of Tq is equal to the height of T. If the sets of outgoing 
edges of so and si are bunches, there is a trivial stable pair (sq, si). If the set 
of outgoing edges of so is a bunch and the set of outgoing edges of si are not 
a bunch (see the left part of Figure [7]), we flip the edge {t,bi,pi) and the red 
edge going out of t. We then call the procedure FlipEdges(^, r) with this 
new red cycle. The root r has now a unique child (si) ancestor of maximal 
state whose set of outgoing edges is a bunch (see the right part of Figure [7j) . 
This call is thus performed at most one time. Finally, if sq is a not a bunch, 
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let (so, bo, go) be a 5-edge with qo ^ r. If qq does not belong to T we flip the 
edge {sq, bo, qo) and the red edge going out of sq. The equivalent automaton 
satisfies Condition (*). It qo belongs to T and is not a descendant of si, 
we flip the edge {t,bi,pi) and the red edge going out of t and we flip the 
edge (so, bo, qo) and the red edge going out of sq. The equivalent automaton 
satisfies Condition (*). If go belongs to T and is a descendant of si, we 
return the automaton together with the edge {so,bo,qo)- 

The previous sequence of flips is described below in the pseudocode 
FlipEdges. The value of previous[r,p] is the pair (so,si), where sq is 
previous[r], i.e. the predecessor of r on its red cycle, and si is the child of r 
which is an ancestor of p. When all states ti are equal to a same state t, the 
code is split into to procedures, FlipEdgesChild and FlipEdgesChil- 
DREN corresponding respectively the case of /? = 1 and the case of p > 1. 

FlipEdges ( automaton A, root r of a maximal tree) 

We assume that the maximal level of the states is positive and that A satisfies Condition 

1 let si, S2, . . , Sp be the children of r which are ancestors of a maximal state 

2 let Pi be a maximal state descendant of (1 < i < p) 

3 let So be the predecessor of r in its red cycle 

4 let r,r2,. -Vk be the ordered list of maximal roots 

5 let T be the tree rooted at r 

6 if there is a predecessor of pi by a 6-edge {U, bi,pi) 

such that U ^ cluster[r] or level[ti] 7^ or (tt G cycle[r] and U ^ [rk ■ ■ r[) 



7 then flip {ti, bi,pi) and the red edge going out of ti 

8 (so) Si) <— previous[r,p] 

9 return A, (so, Sj) 

10 else (ti e[rk..r[) 

11 if there are two distinct states ti (say ti, t2) 

12 then let us assume that t\ < t2 

13 flip (ti, bi,pi) and the red edge going out of t 

14 let T' be the new tree rooted at r 

15 if height[T'] > height[T] 

16 then Sq ^ previous[r] 

17 {so,si) ^ {previous [r],SQ) 

18 return ^, (so, si) 

19 else {height[T'] < height[T]) 

20 flip {t2, b2,P2) and the red edge going out of t' 

21 {s\,S2) previous[r,p2] 

22 return^, (si,S2) 

23 else (all states ti are equal to a state t) 

24 if p > 1 

25 then return FlipEdgesChildren(^, r, t, (sj)) 

26 else return FlipEdgesChild(^, r, i, si) 
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FlipEdgesChild( automaton A, root r of a maximal tree, t, si) 



We assume here that p = 1 

1 let To be the tree rooted at r obtained if we flip {t, 61, pi) and the red edge 
going out of t and keep only r and the subtree rooted at the child sq 

2 if height[To] > height[T] 

3 then flip (t, bi,pi) and the red edge going out of t 

4 Sq ^ previous [r] 

5 {sq,s\) {previous[r\,SQ) 

6 return A, (sq, si) 

7 elseif height[%] < height\T] 

8 then return 

9 else {height[%\ = /iei9/;i[r]) 

10 if the sets of outgoing edges of sq and si are bunches 

11 then return A, (sq, si) 

12 elseif the set of outgoing edges of sq is a bunch 
and the set of outgoing edges of si is not a bunch 

13 then flip (i, 61, pi) and the red edge going out of t 

(we still have height[%] = height[T]) 

14 return FLlPEDGES(^,r) 

15 else (the set of outgoing edges of sq is not a bunch) 

16 let (so, bo, Qo) a 6-edge going out of sq with qo ^ r 

17 if go ^ tree[r] 

18 then flip (sq, bo,qo) and the red edge going out of sq 

19 let / ^ root[qo] 

20 (sq, Si) ^ prewo'us[r',p] 

21 return^, (so,.si) 

22 elseif qo is not a descendant of si 

23 then flip the edge {t, bi,pi) and the red edge going out of t 

24 flip the edge (sq, bo^qo) and the red edge going out of sq 

25 let r' ^ rooi[go] 

26 (so, si) <— premoMs[r',p] 

27 return ^, (so, si) 
28 

29 else [qo is a descendant of si) 

30 return >l, (so, '>o, ^o) 
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FlipEdgesChildren( automaton A, root r of a maximal tree, t, (s^)) 



We assume here that p> 1 and all states ti are equal to a state t 

1 let To be the tree rooted at r obtained if we flip {t, 61, pi) and the red edge 
going out of t and keep only r and the subtree rooted at the child sq 

2 if height[To] > height[T] 

3 then flip (t, bi,pi) and the red edge going out of t 

4 Sq ^ previous [r] 

5 {so,Si) ^ {previous[r],SQ) 

6 return A, {sq, si) 

7 elseif height[%] < height [T] 

8 then flip {t,bi,pi) and the red edge going out of t 

(now height[%] = height[T] and p is decremented) 

9 return FLipEDGES(yl, r) 

10 else {height[To] = height[T]) 

11 if the set of outgoing edges of sq is a bunch and there is an 
integer z > 1 such that the set of outgoing edges of Sj is a bunch 

12 then return A, {sq, Si) 

13 elseif the set of outgoing edges of sq is a bunch 

and the sets of outgoing edges of for z > 1 are not bunches 

14 then flip {t, 61, pi) and the red edge going out of t 

(we still have height[%] = height[T]) 

15 return FLlPEDGES(^,r) 

16 else (the set of outgoing edges of So is not a bunch) 

17 let (so, b, go) a b-edge going out of sq with qo ^ r 

18 if go ^ tree[r] 

19 then flip (so, 60, go) and the red edge going out of sq 

20 elseif go is not a descendant of si 

21 then flip the edge (t, bi,pi) and the red edge going out of t 

22 flip (so, 60, go) and the red edge going out of sq 

23 else flip the edge {t, 62,^2) and the red edge going out of t 

24 flip (so, 60, go) and the red edge going out of sq 

25 let r' ^ root[qo] 

26 (so,Si) ^ previous[r' , p] 

27 return ^, (so, si) 



The function FlipEdges(^, r) returns either an equivalent automaton 
together with a stable pair of states, or an equivalent automaton or an 
equivalent automaton together with a 6-edge. 

After running FlipEdges(^, r) on all maximal roots r, wc get either an 
automaton with a stable pair, or an automaton where each cluster fulfills 
the following conditions 

1. the root of each maximal tree % has a unique child; 

2. there is a 6-edge (xj, bi,yi) such that rj is no more a maximal root after 
flipping this edge and the red edge going out of Xi. 

If the latter case, we flip the b-edge {xi,bi,yi) and the red one starting at 
Xi for all maximal roots but one. We get an equivalent automaton which 
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satisfies Condition (*) as is shown in Figure [HI The pseudocode for this final 
treatment is given in procedure FindStablePair. 

FindStablePair( automaton A) 

1 if the maximal level of states of „4 is 

2 then return FlipEdgesLevelZero(^) 

3 else 

4 A,TZ^ ScanTrees(^) 

5 if i? is a (stable) pair of states (s, t) 

6 then return A, {s, t) 

7 else (i? is empty) 

8 for each maximal root r 

9 do i? <— FLlpEDGES(y4,r) 

10 if i? is a (stable) pair of states {s,t) 

11 then return (s, f) 

12 else {R is a 6-edge) 

13 e[r] ^ R 

At this stage FlipEdges has not returned a stable pair 

14 To ^ a maximal root 

15 for each maximal root r ^ rg 

16 do flip the edge e[r] — {x, b, y) and the red edge going out of x 

17 p <— a maximal node in the tree rooted at 

18 (s,t) ^ previous[rQ,p\ 

19 return A, (s, t) 

5 The complexity issue 

We establish the time and space complexity of our algorithm in the following 
proposition. We denote by k the size of the alphabet A and by n the number 
of states of A. Since A is complete deterministic, it has kn edges. 

Proposition 4. The worst-case complexity of FindColoring applied to 
an n-state aperiodic automaton is 0{kn?) in time and 0{kn) in space. 

Proof. The edges of the automaton can be stored in tables indexed by the 
states and labels. For each state, we compute in linear time two lists of 
predecessor states by red and 6-edges respectively. 

For each state p of the automaton we memorize some data like its cluster, 
its cycle, the root of the tree containing p, the level of p. If level[p\ > 0, 
the data contains also the child of the root of the tree containing p which 
is an ancestor of p. If level[p] = 0, the data contains also the height of the 
tree rooted at p, and the number of p on its cycle. All these data can be 
computed in a linear time with a depth-first search that starts at the states 
that have no predecessors by a red edge. We denote by £ the maximal level 
of the states. 

We first prove that the complexity of FindColoring is at most n times 
the complexity of FindStablePair, when we do not take into account 
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the cost induced by Line [5] of FindColoring for computing the quotient 
automata. Let us show that FindStablePair can be performed in hnear 
time. 

Checking at Line [6] of FlipEdges whether ti € [rfc,..,r[ is done in 
constant time for each ti using the numbers of states at level 0. 

In the worst case, Line [15] of FlipEdges takes a time at most equal to 
the number of states contained in the trees whose roots belong to ]t, ■ ■ ,r[. 
Indeed, in Line [T5l we have to compute the maximal value of height[x] + 
num[r] — num[x] for all x G]t, . . , r[ on the red cycle, where height[x] denotes 
the height of the tree rooted at x. Hence the overall time spent for computing 
FlipEdges(^, r) for all maximal roots r is 0{kn). 

Similarly, in the worst case, FlipEdgesChild(^, r, t, si) takes a time 
at most equal to the number of states contained in trees whose roots belong 
to ]rfc, . . ,r], where is the maximal root before r in the red cycle of r. 
In order to perform the test of Line [2] of FlipEdgesChild, we have to 
compute the maximal value of height[x] + nuni[r] — num[x] for all x €]t, ■ ■ ,r[ 
on the red cycle. In the case where this value is equal to i and the set of 
outgoing edges of sq is a bunch, we flip two outgoing edges of t. We then 
only have to update the data of the states contained in the trees whose roots 
are x ■ ■ ,r] before the second (and last) call to FlipEdges(^, r). It is 
not needed to update num[p] for all states p in the new cycle. It is enough 
to mark the states x € [t, . . ,r[ on the new red cycle (i.e. the states t and 
all states on the path from p to r in the tree of r before the flip). The 
complexity for all calls to FlipEdgesChild is thus at most the number of 
edges, i.e. at most kn. The same result holds for FlipEdgesChildren. 

Moreover, Lines O to [16] of FindStablePair can be performed in time 
at most the number of maximal roots, i.e. at most n. Hence FindStable- 
Pair has a 0(A;n)-time complexity. 

Suppose that the class of a state for the current stable pair congruence 
is stored in an array giving the leader of the class. The Union of two classes 
can be performed in time at most 0(n) while the class of a state is found in 
constant time with the Find operation. 

Since lifting a coloring from an automaton up to another one uses only 
Find operations, it can be done in time 0{kn). 

Finally, we show that the computation of all quotient automata in Line [5] 
of FindColoring during the while loop has an overall quadratic-time com- 
plexity. Since the number of states decreases when two classes are merged, 
the number of merges realized is at most n. Thus the total time for the 
computation of all quotient automata is at most kn?. Note that this com- 
plexity may be improved with a more clever implementation of Union-Find 
operations [_3j. 

Hence FindColoring has a 0(A;n^)-time complexity. The space com- 
plexity is 0{kn). Indeed, only linear additional space is needed to perform 
the computation of the current congruence. □ 
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Figure 2: FlipEdges Lines [6][9l The picture on the left illustrates the case 
corresponding to line [6] of FlipEdges. After the execution of Lines [HE] of 
FlipEdges, i.e. after flipping the edge and the red edge going 

out of ti, we get the automaton on the right. It satisfies the Condition (*), 
i.e. it has a unique maximal tree (here rooted at r). Maximal states are 
colored and the (dashed) 6-edges of the automaton are not all represented. 



6 The case of periodic graphs 

Let the period of an automaton be the gcd of the lengths of the cycles in its 
graph. If the automaton A is an n-state complete deterministic irreducible 
automaton which is not aperiodic, it is not equivalent to a synchronized au- 
tomaton. Nevertheless, the previous algorithm can be modified as follows to 
find an equivalent automaton with the minimal possible rank in a quadratic 
time. 

PeriodicFindColoring( aperiodic automaton A) 

1 .Ao ^ ^ 

2 repeat 

3 A, (s, t) ^ FindStablePair(^) 

4 if there is a stable pair (s, t) 

5 then lift the coloring of A up to 

6 else return 

At Line El it may happen that FindStablePair does not return a stable 
pair of states. In this case, the condition at Line H] is false. At the end of 
PeriodicFindColoring, we get an automaton A that has no stable pair 
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Figure 3: FlipEdges Lines [T2lfT8l The picture on the left illustrates the 
case corresponding to lines [T^ of FlipEdges. The height of the tree T' 
obtained after flipping the edge {ti, bi,pi) and the red edge going out of ti, 
is greater than the maximal level. We get a unique maximal tree rooted at 
r in the same cluster. The picture on the right illustrates the result. 




Figure 4: FlipEdges Lines [T91I221 The picture on the left illustrates the 
case corresponding to the hypothesis of Line [19] of FlipEdges. The height 
of the tree T' obtained after flipping the edge and the red edge 

going out of ti, is not greater than the maximal level. In this case, we also 
flip the edge (i2)&2,P2) and the red edge going out of t2- We get a unique 
maximal tree rooted at r in the same cluster. The picture on the right gives 
the resulting cluster. 
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Figure 5: FlipEdgesChildren Lines [T^lfT^ The picture on the 
left illustrates the case corresponding to the hypothesis of Line [13] of 
FlipEdgesChildren. Let % be the tree rooted at r obtained if we flip 
{t,bi,pi) and the red edge going out of t and keep only r and the subtree 
rooted at the child sq. The states of the tree Tq rooted at r are represented in 
salmon in the left part of the figure. The state sq is a bunch. After hipping 
the edge {t,bi,pi) and the red edge going out of t, we get the automaton 
pictured in the right part of the figure. The The tree Tq (new tree Tq) is 
the tree rooted at r obtained if we fiip {t,b2,P2) and the red edge going out 
of t and keep only r and the subtree rooted at the child si. Its states are 
colored in salmon. Its height is the same as before. 
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Figure 6: FlipEdgesChildren Lines [2511771 The picture on the 
left illustrates the case corresponding to the hypothesis of Line [23] of 
FlipEdgesChildren. Let % be the tree rooted at r obtained if we flip 
{t,bi,pi) and the red edge going out of t and keep only r and the subtree 
rooted at the child sq. The states of the tree % rooted at r are represented 
in salmon in the left part of the figure. The state sq is not a bunch: it has 
a 6-edge {so,bo,qo) with qq = si. After flipping the edge (i, 62,^2) and the 
red edge going out of t, and flipping (so,6oi'?o) and the red edge going out 
of So, we get a unique maximal tree rooted at r in the same cluster (see the 
right part of the figure) . 
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Figure 7: FlipEdgesChild Lines [T2]fTil The picture on the left illustrates 
the case corresponding to Line[T2]of FlipEdgesChild. After the execution 
of lines [T2]fTil of FlipEdgesChild, i.e. after flipping the edge (t, and 
the red edge going out of t, we get the automaton on the right. The root r 
has a new single child si ancestor of a maximal state, whose set of outgoing 
edges is a bunch. The new tree rooted at r has here the same level ^ = 2 as 
before and FlipEdges(^, r) is called a second and last time. 




Figure 8: FindStablePair Lines [TUfTSl The picture on the left illustrates 
the case corresponding to the conditions at Line [U] of FindStablePair. 
After the execution of FindStablePair, two edges of the predecessors in 
the red cycle of all but one maximal roots are flipped. The new cluster is 
pictured on the right part of the figure. It has a unique maximal tree. By 
Lemma E] the pair (6, 15) is stable. 
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(it is made of a cycle where the set of outgoing edges of any state is a bunch) . 
Lifting up this coloring of to the initial automaton Aq leads to a coloring of 
the initial automaton whose minimal rank is equal to its period. 

This result can be stated as the following theorem, which extends the 
road coloring theorem to the case of periodic graphs. 

Theorem 5. Any irreducible automaton A is equivalent to a an automaton 
whose minimal rank is the period of A. 

Proof. We prove the result by recurrence on the number of states of A. 

Let us assume that A is equivalent to an automaton A' which has a 
stable pair (s,t). Let B' be the quotient of A' by the congruence generated 
by Let d be the period of A' (equal to the period of A) and d' the 

period of B'. It is clear that d' divides d. Let w be the label of a cycle 
around a state s oi B' . Let {pi, . . ,Pk} be the congruence class of the states 
of A! that corresponds to s. Hence, after renumbering the states of s, there 
is an integer i and a cycle in A! of length \w\ x i: 

www 
Pi ^P2 ^ ■■ ^ Pi+1 = Pi- 

Let us assume that i > 1. Since {pi,P2) is a stable pair, there is a word u 
such that pi ■ u = p2 ■ u. Since A' is irreducible, there is a word v labeling a 
path from pi ■ ii to pi. Thus there are cycles around pi with labels wuv and 
uv. This implies that d' divides \w\ and thus that d = d'. This still holds if 
i = 1. 

Suppose that B' has a coloring B" which has rank d. We lift this coloring 
up to an automaton A" equivalent to A" . Let us show that A" has rank d. 
Let / be a minimal image of A" and J be the set of classes of the states of 
/ in B". The set J is a minimal image of B" . Two states of I cannot belong 
to the same class since / would not be minimal otherwise. As a consequence 
/ has the same cardinal as J. 

Let us now assume that A has no equivalent automaton with a stable 
pair. In this case, we know that A is made of one red cycle where the set of 
edges going out of any state is a bunch. Since A has rank d, the length of 
the cycles is d. Hence A has minimal rank d. □ 

Since the modification of FindColoring into PeriodicFindColoring 
does not change its complexity, we obtain the following corollary. 

Corollary 6. Procedure PeriodicFindColoring finds a coloring of min- 
imal rank for an n-state irreducible automaton in 0{kn'^). 

Acknowledgments The authors would like to thank Florian Sikora and Avra- 
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